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Abstract 



In this paper, we consider the problem of model equivalence for quantum systems. 
Two models are said to be (input-output) equivalent if they give the same output for 
every admissible input. In the case of quantum systems, the output is the expectation 
value of a given observable or, more in general, a probability distribution for the result 
■ of a quantum measurement. We link the input-output equivalence of two models to 

' the existence of a homomorphism of the underlying Lie algebra. In several 

, Cartan decomposition of the Lie algebra su{n) is useful to find such a homomorphism 

04 ■ and to determine the classes of equivalent models. We consider in detail the important 

' cases of two level systems with a Cartan structure and of spin networks. In the latter 

, case, complete results are given generalizing previous results to the case of networks 

of spin particles with any value of the spin. In treating this problem, we prove some 
I instrumental results on the subalgebras of su{n) which are of independent interest. 
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Dynamical models of quantum systems have been recently the subject of investigation, 
concerning their structural properties, by use of methods of control theory. Appropriate 
definitions of controllability and observability of quantum systems have been given and 
practical conditions to check these properties have been proposed (see e.g. P, [S], [TU]. 
' E])- In many cases, the tools used are the ones of Lie algebra and Lie group theory. 

^ . Information on the properties of the dynamics is obtained by a study of the structure of 

a Lie algebra associated to the system and how this relates to the particular equations 
at hand. This geometric approach has proved useful not only to analyze the dynamics 
but also to design control laws. This approach can also be used to study problems of 
parameter identification of quantum systems and this is the subject of the present paper. 
' In particular, the problem we shall study is the classification of models of quantum systems 

whose behavior cannot be distinguished by an external observer. We shall call these models 
(input- output) equivalent. This problem is motivated by several experimental scenarios. In 
particular consider a molecule which is a network of particles with spin with all the other 
degrees of freedom neglected. A model Hamiltonian is associated to this system in which 
parameters modeling the interaction between particles as well as the interaction with an 
external electro-magnetic field are unknown. Also, the initial state of the system might 
be unknown. In experimental scenarios such as Nuclear Magnetic Resonance and Electron 
Paramagnetic Resonance, it is possible to drive the system with a magnetic field and measure 
the expectation value of a given observable, for example the total spin in a given direction. 
The question of fundamental and practical importance is to what extent, with this type of 
experiments, it is possible to distinguish between different models. As we shall see in this 
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paper, this question is related to the existence of a particular Lie algebra honiomorphism 
which relates the equations of the two models. 

The main results of this paper are the solution of the model equivalence problem for a 
class of two level systems in Theorem|5J a number of auxiliary results (Theorems I4I6|I on the 
structure of the Lie algebra su{n) and its subalgebras and Theorem [7| where we completely 
solve the problem of characterizing equivalent models for networks of spin. The latter result 
generalizes results previously obtained in [2] and [7] , which were proven only for networks 
of spin i and I's, to networks of interacting spins of any value and where the spin itself is 
an unknown parameter to be identified. A further motivation for this research can be found 
in J5 where it was shown that thermodynamic methods commonly used to identify the 
parameters of spin networks such as in molecular magnets 0, [H] are not always adequate. 
The generalization presented in this paper is obtained through a Cartan decomposition 
technique recently presented in j2| which helps determining the homomorphism between 
equivalent models in the form of a Cartan involution. 

The paper is organized as follows. In Section 13 we describe the problem of model 
equivalence for quantum systems. In Section |31 we link the equivalence of two models to 
the existence of an appropriate Lie algebra homomorphism. In several cases the structure of 
the dynamics is related to a Cartan decomposition of su{n) and suggests the form of such a 
homomorphism as well as of the classes of equivalent models. We give a two level example 
in Section 01 and treat the case of general spin networks in Section [3 Instrumental to the 
solution of the model equivalence problem for spin networks are some results of independent 
interest concerning the existence of subalgebras of su(n) with specific features. The proofs 
of these results are presented in Section |B| Concluding remarks are given in Section |7| 



2 The problem of model equivalence for quantum sys- 
tems 

Consider a model Hamiltonian for a quantum system, H{t) :— H{u{t)), where, in a semi- 
classical description, the dependence on time is due to the interaction with external fields, 
u := u{t), which play the role of controls. The evolution of the state of the system, described 
by a density matrix p := p{t), is determined, other than by H, by the initial state p(0) — po- 
In particular, p is the solution of the Liouville's equation 

P=HH,p], (1) 

with initial condition p{0) = pQ. For this system, we measure an observable S. Considering, 
for simplicity, a Von Neumann-Liiders measurement'^, writing S in terms of orthogonal 
projections 

S:=J2X,U„ (2) 

j 

the probability of having a result Xj, when the state is p, is given by 

Pj ■■= Tr{U,p). (3) 

As the probabilities Pj are the only information that can be gathered by an external observer, 
it is motivated to ask what classes of models {H{u),po} will give the same probabilities, 
for any functional form of the control it. In other terms, we ask what classes of models are 

^Natural extensions can be made to general measurements |4| for the related issue of observability of 
quantum systems. 



2 



indistinguishable by experiments that involve driving the system with controls, in a given set 
of functions, and measuring a given observable. These models will be called (Input- Output) 
Equivalent. 

It is appropriate to treat the case where the expectation value of the measurable S, i.e. 
the 'output', 

y ■■= Tr{Sp), (4) 

is measured. Not only this is the case in several experimental situations, such as nuclear 
magnetic resonance, but it is not a significant restriction as compared to the case where 
the probabilities Q are considered. As the structure of the output I^J is the same as the 
one of the outputs ||2J), the passage from the treatment for the expectation value to the one 
for probabilities corresponds to extending a single output treatment to a multiple output 
treatment. This can be accomplished without difficulties. 

In order to render the problem of characterizing the classes of equivalent models treatable, 
we need to assume some structure on the Hamiltonian H . This corresponds to the passage 
from unstructured uncertainty to parametric uncertainty often discussed in identification 
theory (see e.g. In particular, it is often the case that the Hamiltonian H = H{u) has 

the bilinear form 

m 

H := Ho+Y,HjUj{t), (5) 

for some control functions ui, ...,Um, and internal Hamiltonian Hq and interaction Hamil- 
tonians Hj^s, j — 1, ...,m. These can be considered as Hermitian matrices of dimension n, 
i.e. elements of isu{n), where n is the dimension of the system, assumed finite. Moreover, in 
many cases, Hq and the iJj's belong to two orthogonal complementary subspaces of isu{n) 
corresponding to a Cartan decomposition of su{n) |^. These are two subspaces iK. and iV, 
such that the subspaces of su{n), K, and V, satisfy the commutation relations 

[/c,/c]c/c, [k:,v]<zv, [p,v]cic. (6) 

If the system is a multipartite system, every Hj is a linear combination of Hamiltonians 
modeling the interaction of each individual system with the external field. In matrix nota- 
tion, Hj is a linear combination of elements of the type lig)!®---®!®^®!®---®!, where 
L is an Hermitian matrix of appropriate dimensions and all the other places are occupied by 
identities 1. Also, Hq is very often a linear combination of elements modeling the interaction 
between two subsystems, which can be written as tensor products of matrices equal to the 
identity except in two locations. In these cases the relevant Cartan decomposition for the 
Hamiltonian |(SJ) can often be chosen of the odd- even type described in the recent paper [5]. 
Also if 5 is a sum of observables on each individual subsystem, i.e. total angular momentum 
(see e.g. 52 )i it can always be written as sum of tensor products all equal to the identity 
except in one position. In these cases iS, belongs to one of the subspaces of the Cartan 
decomposition^ . 

In the following, we shall consider, as standing assumption, only finite dimensionality of 
the Hamiltonian H and the bilinear form Q and will make precise the assumptions on the 
Cartan structure of the Hamiltonian when needed. 

* Notice that the situation may be different if we consider the case of a single output given by the 
expectation value J3} and the case of several outputs given by the probabilities in • 
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3 Model equivalence and Lie algebra homomorphisms 

Consider two models with Hamiltonian of the form jSJ and output of the form Q 

m 

P = [-ii.HQ + ^ HjUj), p], 

m 

The foUowing theorem hnks the existence of an appropriate Lie algebra homomorphism to 
the equivalence of the two models. 

Theorem 1 Let n and n' be the dimensions of the two models respectively. Let cj) 

he a homomorphism, (j) : u{n) — > u(n'), and (j)* its dual with respect to the standard inner 
product < A,B >:— tr{AB*). Assume 

-iH[,^(l){-^iHo), -iH'j = (j){-iH.j), (j)*{iS')^iS. (9) 

Then ifipQ — (t>{ipo) the models are equivalent. Viceversa, if the models are equivalent and 

0) is observable, then 

ip'o = 0(zpo). (10) 

Proof. Multiply {T)) and ^ hy i and then apply (j) to the equation obtained from ([7|). 
Combining the two resulting equations, we obtain 

j^iip' - (/.(zp)) = m~iHo) + ^ ^{-tH,)u„ip' - (11) 

j 

Now, if (|10|l is verified, then ip'{t) — (f){ip{t)), for every t and for every control. Therefore 
we have 

Tr{S'p') = Tr{-iS'i{p')) = Tr{~iS' (t>{ip)) = Tr{4>*{~iS')ip) = Tr{Sp), (12) 

and the two models are equivalent. Viceversa, assume the two models are equivalent. From 

1) 12|l . we have 

TritS'iip' -<t>{tp))it))^0, (13) 

for every t. Writing the solution of |[TT|) as {ip' - <j){ip){t) = X{ip' - (j){ip)){0)X* , where X 
is the solution of the (Schrodinger) operator equation X = {(f){—iHo) + J^j 4>{—iHj)uj)X, 
X{Q) = 1, we have 

Tr{X*iS'X{ip'^ ~ 0(^po))) = 0. (14) 
As the system ((H)) is observable, we have that X*iS'X span all of su{n'), which implies 

As we shall show in the remainder of the paper (cf. also PJ), it is possible for cases of 
physical interest to give a stronger version of Theorem^ In particular, it is possible to show 
that the existence of a homomorphism satisfying (|5J| is also necessary for equivalence of 
two models. This way, we can characterize all the classes of equivalent models in terms of 
homomorphisms. We shall do this for a two level example in the next section and general 
spin networks in Section |S1 In both cases we exploit a Cartan decomposition underlying the 
dynamics of the models. In general, more structure will have to be assumed to avoid trivial 
cases. For example, if S* = S" is a scalar matrix, then every two models are equivalent. To 
avoid this case, a reasonable extra assumption is the observability of the two models. Also, 



p(0)=po, y^Tr{Sp), (7) 
p'{Q)^p',, y' = Tr{S'p'). (8) 
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we need to assume that the initial states are not both perfect mixtures, otherwise, with 
S — S' , the output for any two equivalent models will be the same, independently of the 
dynamics. Moreover —iHj and —iHj, j = 0, . . . , m, may be in general assumed traceless, 
as the trace only adds an extra common phase factor to the dynamics, which cannot be 
detected. 



4 Model equivalence of two level systems 

Consider a spin i particle which is driven by an electro-magnetic control field along the 
z axis, interacts with a constant unknown magnetic field along a (unknown) direction in 
the X — y plane and has unknown initial state. The practical question is to what extent, 
by driving the system with the control field and measuring the average value of the spin 
magnetization in the z direction, it is possible to obtain information about the unknown 
parameters of the system. This type of model has a Cartan structure which is shared by 
several other models of physical interest and is instrumental in finding a homomorphism 
between equivalent models. We describe this below. 

The Lie algebra su(2), which is the relevant Lie algebra in the two level case, has, up to 
conjugacy, only one Cartan decomposition (or in other terms all the types of decompositions 
coincide in this case) which corresponds to the classical Euler decomposition of the Lie group 
SU{2) 0. This extends to a decomposition of u{2) which can always be written, as 

u{2) = IC®V. (15) 

Here K- and V satisfy the commutation relations in ((HJ and are given, up to conjugacy, by 

K. := span{i(Tz}, V := span{iax,icry,il2x2}, (16) 

where l2x2 is the 2x2 identity matrix and ax, Uy and Uz are the Pauli matrices 

1 /O l\ I ( Q i\ 1 /I \ 

^-^=2(,1 oj' ''y^■^2\-^ oj' ^^^^=2 1,0 -1 j ' ^^^^ 

The dynamical and output equation, for the above model of a spin ^ particle in an electro- 
magnetic field, can be written as 

p=[A + ia:,u{t),p], y^Tr{(Tzp), p(0) = po, (18) 

where po is an unknown initial density matrix and A :— xiux +yiay, with x and y unknown. 
This model has a Cartan structure in that is in T' and iaz (the control and observation part) 
is in /C, with /C and V defined in p6|l . We assume x^ + y^ which implies controllability 
and therefore observability for this model. The following result characterizes all the 
classes of equivalent models in terms of Lie algebra homomorphisms. 

Theorem 2 Consider two models 

p=[A + iazu{t),p], y^Tr{(Tzp), p{0) = po, (19) 

p' = [A' + iaMt),p], y = Tr{azp'), p'{0) = p'o, (20) 

with Po and p'q not both equal to scalar matrices (representing perfect mixtures) and A and 
A' given by 

A :— xiax + yic^y, and A' :— x'iux + y'icTy. (21) 
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Assume 

+ ?y V and x'^ + y'^ / 0. (22) 

Then the two models are equivalent if and only if there exists an automorphism (j) : u{2) — > 
u(2) with 

^*{ia,)^ia, (23) 

and 

A' = (j){A), (j){ia,) ^i(7„ ip'Q^4>{ipo). (24) 

Proof. It is clear that if the automorphism exists, satisfying (|23() H24() . the two models 
are equivalent. This follows from a direct application of Theorem^ To prove the opposite, 
first notice that, from the equivalence assumption, we have 

y{t) := Tr{a,p{t)) = Tr{a,p'{t)) := y\t), (25) 

for every t > and every admissible control. 

We consider an inner automorphism (j) of the following type 

(j}{L):^e-'°"'^Le'°"'% Leu{n), (26) 

as a varies in R . 

Clearly H23|) and the second one of H24I) are verified for any a G R . Moreover 



with 



and 



4i{A) = xia^ + yiay, (27) 



'''' 



:= I ""i"! I ^ (29) 



Also, if we write 



we have 



-sin{a) cos{a) 

ipo Pxi<Jx + PyicTy + PziCTz + 

ipQ := p'Ja.j, + p'yicTy + p'^ia^ + \il 



(30) 



(ipo) = Pxil^x + PyiCy + Pzi^z + ^il, (31) 



' > - f ) . (32) 



with 

\Py J \Py 

Using the equivalence assumption H25|l at t = we immediately obtain 

Pz = p'z- (33) 

Moreover, differentiating 1)25(1 using the dynamical equations ((19(1 ((20(1 . we obtain 

Tr{p[a,,A])^Tr{p'[a,,A']). (34) 

Writing this at time t — and using the definitions ((21(1 and ((30(1 along with the commutation 
relation for the Pauli matrices 

[iaj;,iay] = iaz, [iay,iaz] ^ icr^, [iaz^iox] = iay, (35) 
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we obtain 



PyX - pxy = p'yx' - p'^y'. (36) 

Differentiating H34|l and using the fact that the resulting equation has to be vaHd for every 
value of the control, we obtain the two equations 

Tr{aM, [A,p]]) = Tr(a,[A', [A',p']]), (37) 

and 

Tr{ia,[A, [a,,p\]) = Tr{ia,[A' , [a„p']]). (38) 
From equation H38|l. as for equation (|36|l . we obtain 

xp^ + ypy^x'p'^+y'p'y. (39) 
From equation H37|l. we obtain 

{x^ + y^)Tr{a,p) = {x'^ + y'^)Tr{a,p'). (40) 

Using the fact that Tr{azp) is not always zero (because of the controllability condition l^'22\i ) 
^ and equation 125|l . we have 

Therefore, for some a, we can write 

with Ka in (|29|l . and this, compared with (|28|l and H27|l . gives the first one of (|24|l . To obtain 
the third one (with the same 0), we recall (I22II that x'^ + y'^ ^ 0. Letting (^^^ 
and using (|43() . we can write H36|l and (|39|l . respectively, as 

[a^, y]J[p.,Pyf = [x, y]Klj[p',,p'f, (44) 

[xM[p.,PyV = [^.y]Kl[p'^,p'yV. (45) 
Since iCj^ commutes with J, we can write these as 



Since + = —det | ' | ^ 0, we can write 



= (47) 



^From controllability I22i . we cannot have 

Tr(a^pit)) = Tr{(7^p'{t)) = 0, (41) 

for every control. This would mean that, for every reachable evolution operator X, solution of the 
(Schrodinger) operator equation X = [A + iazu)X, X'a^X would be orthogonal to pp. However, be- 
cause of controllability X may attain all the values in SU{2) and therefore X*crzX span, as X varies, all of 
isu{2). Therefore, X'tr^X is always orthogonal to po only if po is a multiple of the identity, which we have 
excluded. 
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and therefore 

[p'x.p'v]^[P-^Pv]^ (48) 

which along with p'^ — pz gives 

ip' = cbiip). (49) 
This condudes the proof of the Theorem. □ 

5 Model equivalence of spin networks 
5.1 Set of models of spin networks 

We consider a network of n particles with spin that interact according to Heisenberg inter- 
action. In particular, we denote the spin of the j*'* particle by Ij and by Nj :— 2lj + 1 the 
dimension of the Hilbert space for the state of the j*'* particle. The dimension of the Hilbert 
state space associated to the entire network is TV :— Y[j=i -^r '^^^ class of Hamiltonians we 
consider are of the form 

H{t) -.^ i{A + B^u^it) + ByUyit) + Bzu,{t)), (50) 

where A, modeling the Heisenberg interaction among the particles, and Bx,y,z, modeling the 
interaction with external fields, are given by 

A := —iJ2k<l.k.l=l Jklihx.lx + Ikydy + Ikz,lz)- 

(51) 

Bv ■= -i{YJk=ilkIkv), for v^x.y, or z, 

respectively. Here and in the following we denote by Ikivi,...,kr.vr^ for 1 ^ ki < ■ ■ ■ < kr < n 
and Vj 6 {x, y, z}, the N x N matrix which is the Kronecker product of n matrices where in 
the j*^ position we have the Nj x Nj identity if j ^ {ki, . . . , kr}, while if j = kg we have the 
Nj X Nj representation of the Vs component of spin angular momentum for a particle with 
spin Ij. Such matrices are given by the Pauli matrices 1)17(1 in the case where Ij = ^ and 
can be calculated for every value of the spin (see e.g. J2] Section 3.5). With some abuse 
of notation, we shall continue denoting these matrices by ax, <Jy and cr^, without explicit 
reference to the value of the spin. These matrices have several properties we shall use in the 
following. In particular, they satisfy the commutation relations (|35|l and 

al+Gl+al= + 1)1n,xn, (52) 

(see e.g. formula (3.5.34a) in [T^). 

The real scalar parameter Jki in (|51|l is the exchange constant between particle k and 
particle I and the real scalar parameter jk is the gyromagnetic ratio of particle k. We assume 
that the spins of the network have all non-zero and different gyromagnetic ratios. We can 
associate a graph to the model, where each node represents a particle and an edge connects 
two nodes if and only if the corresponding exchange constant is different from zero. It is not 
difficult to see that if the model is controllable then, necessarily, this graph is connected. 
Moreover, controllability implies observability for every output of the form where S* is a 
non scalar matrix . In our case, we assume to measure the expectation values of the total 
magnetization in the x, y, and z direction, given as in where S is one of the matrices: 

n 

Sy := ^/fet,, with V e {x, y, z}. (53) 
fc=i 
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A model of the type above described will be denoted by S := S(n, Ij, Jki,lk, Po)j where the 
parameters n,lj, Jkhjk, po, which determine the model, are unknown. We will assume to 
have two controllable models E and S' :— Z^-, J(,;, 7^., pp) which satisfy the previous 

requirements and we look for necessary and sufhcient conditions for these two models to be 
equivalent. We shall mark with a prime, ', all the quantities concerning the system S'. 

5.2 Relevant homomorphism of su{n) 

In 8 a method was described to construct a Cartan decomposition of the Lie algebra su{N) 
for a multipartite system, starting from decompositions of the Lie algebras su{Nj) associated 
to the single subsystems, each of dimension Nj, with N := 11^=1 -^j- particular, we have 
the following result. 

Theorem 3 fjl^. Section 5) Consider a multi partite system with n subsystems of dimen- 
sions iVi, . . . , Nn- Consider the Lie algebra u{Nj) related to the j-th subsystem and a Cartan 
decomposition 

u{Nj)^JCj®Vj, (54) 

of the type AI or All ^ . Denote by aj (Sj) a generic element of an orthogonal basis of ilCj 
(iVj). Let the (total) Lie algebra u{NiN2 ■ ■ ■ Nn) be decomposed as 

iu{NiN2---Nn)^Io®Ie- (55) 

2o (Xe) is the vector space spanned by matrices which are the tensor products of an odd 
(even) number of elements of the type CFj. Then u{NiN2 ■ ■ ■ Nn) = ilo ffi is a Cartan 
decomposition i.e. 

[iIo,iIo\ C ilo, [Ho, He] C He, [He, He] ^ Ho. (56) 

The decomposition H55II is called a decomposition of the odd- even type. 

Associated to a Cartan decomposition H56|l is a Cartan involution <p which is the identity 
on Ho and multiplication by —1 on He- A Cartan involution is clearly a homomorphism. 
The structure of system (|5()|l and H51(l suggests that it is possible to choose this Cartan 
involution as a homomorphism mapping the equations of two equivalent models as in (|5J). 
In fact, assume that there is the same number of subsystems (spin particles) in the two 
models and that corresponding subsystems have the same dimension (namely the same 
spin). If we can display a decomposition 1)54(1 of the type AI or All for every (spin) su{Nj), 
such that i(Jx,y,z G ^j, then, for every value of the parameters, it holds that B^^y^zi') G Ho 
and A{') G He. As shown in the following Theorems I4I6I decompositions of this type exist. 
We shall see in the following subsection that the Cartan involution associated to an odd-even 
type Cartan decomposition is the correct homomorphism to describe classes of equivalent 
spin networks. In fact, not only models which are related by such a homomorphism are 
equivalent (according to Theorem^ but the opposite is true as well. In other terms, two 
equivalent models are either exactly the same or are related through such a homomorphism. 

The following three Theorems show the existence of a decomposition of su{Nj) of the type 
AI or All where the subalgebra ICj contains the matrices ia^, i(Jy and icr^. Equivalently, 
they show the existence of a subalgebra of sp{^) (type All) or so{Nj) (type AI) conjugate 
to the Lie algebra spanned by iox, ioy and ia^. The proofs are presented in the following 
section. We shall see that the situation is different for integer and half integer spins. 

®In a decomposition AI ICj = so{Nj) and Vj = {so{Nj))-^ up to conjugacy. In a decomposition All 
Kj = sp{Nj/2) and Vj = {sp(Nj /2))-^ up to conjugacy. 
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Theorem 4 // the dimension Nj of the system is even (half integer spin (Fermions)) there 
exists a suhalgehra of sp{^) conjugate to the Lie algebra spanned by iax, icTy and iaz- 

Theorem 5 // the dimension Nj of the system is odd (integer spin (Bosons)) there exists 
a subalgebra of so{Nj) conjugate to the Lie algebra spanned by ia^, i^y cind iaz- 

Theorem 6 // the dimension Nj of the system is even (half integer spin (Fermions)) there 
is no subalgebra of so{Nj) conjugate to the Lie algebra spanned by iax, io'y and iaz- 

5.3 Necessary and sufRcient conditions for model equivalence 

In this subsection we will prove the equivalence result concerning models of spin networks. 
This is given by the following Theorem. 

Theorem 7 Let E :~ "E^n^lj, Jki,^k, Po) o,nd S' :— S](n', Z^-, J^^, 7^, Pq) be two fixed models 
(see equations il5fl)) . ). Assume that both models are controllable, that for model S (^'), 
all the "fk ("f'k) are non-zero and different from each other, and that po and p'q are not both 
scalar matrices. Then E is equivalent to E' i.e.: 

y,{t):=Tr{S,p(t))=y'St):=Tr{Sy{t)), for v(,{x,y,z}, (57) 

and for every control Ux,Uy,Uz, if and only if the following condition holds: 

Cond. (*).• 

L n — n' 

Up to a permutation of the of the set {!,..., n} (i.e. a permutation of the 
indices for the particles) 

2- Ik = ik, 
and 

3. 

h = I'k- (58) 
4- One of the following two conditions holds 
(a) 

A^ A', and po = p'q (59) 

(b) Given the Cartan involution cj) associated to the decomposition of the 
odd-even type as from Theorem\^ 

A' = 4){A), and ipj, = (f){ipo) (60) 

The Theorem says that, under appropriate controllability assumptions, two equivalent 
models for spin networks are equivalent if and only they have the same number of particles, 
corresponding particles have the same spin, and their dynamical model and initial state are 
either the exactly the same or are related through the Cartan involution associated to a 
decomposition of the odd-even type. In practical terms, given a general spin network, by 
driving the network with an external electro-magnetic field and measuring the total spin in 
the x, y and z direction, it is, in principle, possible to identify the number of particles, their 
spin, the gyromagnetic ratios of every spin and the exchange constants only up to a common 
sign factor, if the initial state is not known. The proof that Condition (*) implies equivalence 
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is an application of the general property of Theorem^ The proof that equivalence implies 
Condition (*) is considerably longer. However, several results can be obtained with proofs 
that are formal modifications of the ones presented in j2| for the special case of spin ^ 
particles. We shall focus on the new part of the proof needed to generalize to the case of 
unknown spins. 

Condition (*) implies equivalence 

It is clear that if holds, then the two models differ possibly only by a permutation 
of the indices of the particles. So they are equivalent. Assume now that Condition (*) holds 
with (|60|l and assume for simplicity (and without loss of generality) that the permutation 
of indices is the identity. Let be the Cartan involution associated to the decomposition of 
the odd-even type. We notice that 

(l)*{iSy) ^ iSy ^ iS'^, v = x,y,z. (61) 

In fact, given any C G u{N), we can write C — Co + Ce, with Co G Ho and Ce £ He, it 
holds: 

Tr(0*(i5„)C) ^r((^5„)0(C)) = Tr{{tS^){Co - Ce)) = Tr{{tS^)Co) = Tr{{iS,)C), 

which, since it has to hold for every C, gives (I61|l . Equations (|60|l and H61() imply that 
equation @ of Theorem ^ holds. Since we also have ifTUIl . from (|^ . we conclude that the 
two models are equivalent using Theorem ^ 

Equivalence implies Condition (*) 

The technique used in |2| to prove this result for network of spin 1/2 particles extends to 
the general case treated here. However further analysis is required in this case, in particular 
to prove that equivalent spin networks have the same values of the spins, while in |2] it was 
assumed that the networks were composed by all spin i's. The main reason why the proof 
in can be extended to this case is that the basic commutation relations, which were the 
essential ingredient of the proofs in [2] still hold. More precisely, the matrices ax , cTy and CTz 
still satisfy, for every value of the spin, the commutation relations H85|l . This fact implies 
that it also holds: 

if k ^ {ki, . . . , kr} 
if 3j with k = kj and v^. — vj , 
if 3j with k — kj and v^. ^ Vj 

(62) 

Assume now that the two models E and S' are equivalent. Then, using exactly the same 
arguments as in the proof of Proposition 4.1 of P], we obtain that the number of the spin 
particle must be the same, namely n = n' , and, up to a permutation of the indices, 7^ = 7^., 
Vfc G {1, which is part 1. and 2. of Condition (*). Moreover as in Proposition 4.1 of 

P], we obtain 

Tr{h.p{t)) ^Tr{I',^p'{t)),yk e {l,...,n},yv e {x,y,z}. (63) 

Here I'f,^ is defined as Ikv but for S' and, at this point, it may be different from I^v since 
we have not shown yet that corresponding spins must be equal. To prove this fact, we shall 
use Lemma |5 . II below. The proof of this Lemma is a generalization of the proof of Lemma 
5.2 in [2j where we use the general property 152(1 instead of the corresponding property for 
spin i's. We postpone this proof to Appendix A. 



'felDl,...,fcrllr J J-kvr. 






kivi^...,kj [vjV^] ,...,krVj^ 



independently of the values of the spins. 
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Lemma 5.1 Assume that for all t > 0, all possible trajectories p{t) of S and corresponding 
p'{t) of E', for fixed values 1 < ki, . . . , kr < n, vj ^ {x, y, z} and for given constants /3 and 
/?', we have: 

PTr {h,.,,...,Kv.p{t)) P'Tt (4,,,...,fc,„,p'(t)) . (64) 

Then 

1. For any pair of indices k,d d {1, . . . , n} with k G {fci, kr} and d ^ {fci, fc^}, 

for any value v G {x, y, z}. 

2. For any pair of indices k,d both in {fci, kr}, (for example k — ki, d = ^2) then 

(66) 

In words, formula Ht)6|) means that from (|64|l . it is possible to derive a new formula as 
follows. Select two indices in the set {fci, fcr}, k and d. One of the two indices (say d) 
disappears from the subscript in the matrices / and corresponding /'. However a coefficient 
^di^d + 1) s-iid l^Q-d + 1) appears in the left and right hand side, respectively, as well as a 
coefficient Jj,j and J^^. 

We shall now prove that, under the assumption of equivalence, the squares of the ex- 
change constants Jdk and J^j, must be proportional, with a proportionality factor common to 
all pairs of indices d and fc and this will also be instrumental in the proof of 3. of Condition 
(*)• ^ 

Fix any 1 < fci < fc2 < n, then, by applying statement 1. of Lemma |5.1I i.e. equation 
l)65|l with fc = fci , d = fc2 to equation (|63() with fc = fci , we have: 

Jk,k2Tr{h,vuk2V2P{t)) = J'k^k^Tr [l'k,v,,k2V2P' i^)) ' ^ ""I'^a e {x,y,z}. (67) 

Now, to the previous equality, we apply statement 2. of Lemma l5. II i.e. equation ]<j<i\i with 
fc = fci and d = fc2 to get: 

{Ik2{lk2 + l))Jlk2Tr{Ik,.,p{t)) - (4,(4, + l))4?fc,rr (4,,p'(0) , 

which, by equation H()3|) . implies: 

(/fc,(?fc, + 1))J|^,, - (4,(4, + l))jZk2- (68) 

Using the facts that the two indices fci and fc2 above are arbitrary and that the graph 
associated to the network is connected, by the controllability assumption, it is easy to see 
that there exists a positive constant a G R such that, for all 1 < d < fc < n: 

Jlk^a'j'lk and + 1) = + 1). (69) 

Using 169|) . we can now prove 3. of Condition (*). We will do this using some lemmas 
and arguing by contradiction. First notice that from (|69|) . we have that if there exists a 
fc G {1, . . . , 71} such that Ik = 4.1 then necessarily — 1, thus Ij = 4 for all j = 1, . . . , n, 
namely all the particles have the same spin. So if we assume that (|58|l does not hold, without 
loss of generality, we can assume Zi > l[. Using equation H69|l . we get that Ij > I'j for all 

j = 1, . . . , n, thus also Nj > N'-. Let i? := ^ = nj=2 ^3 and i?' := fj- = n"=2 
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Lemma 5.2 For alH e R , and all the admissible trajectories p and corresponding trajec- 
tory p', we have: 

Tr ((e-^* lii.xii')4(e-*"^* ycR')p' [s)) , ^^^^ 

for aU s > 0. 

Proof. First we notice that from the Campbell-Baker-Hausdorff formula, we have: 

fc=0 

and an analogous equation for E'. Moreover, by applying Lemma l7. II in Appendix A, with 
W — Iiv, W — and fc = 1, V = z, we have: 

Tr{ad,^^g,i^^^Iiyp{s)) = Tr {ad,„^^i^,^^,I[^p'{s)j . 

Now we can apply again Lemma l7.1l to the previous equality, to get: 

Tr (ac'L®ifixfl^i"'^(*)) = Tr {a(f^a,mn'xR'^'ivP'^^^) ' 
By applying repeatedly this procedure we obtain: 

Tr (H>,®ifixfl^i"^(*)) = Tr {adt^^i^,^^,l'ivP'{s)^ , 

for all fc > 0. Using this in (|71|l . equation H70|l follows. □ 

The proof of the following lemma is given in Appendix A 
Lemma 5.3 The following formula holds: 

(e*"^* ® lflxfl)/ix(e-''"^* ® Iflxi?,) := PjVi(i) <E> Irxr, (72) 

where the matrix Pni{') is periodic with period 2tt. Moreover 

PN,{n) = -Pn^{0) = -<J^. (73) 

Using Lemmas 15 . 21 and 15 . 31 we can now conclude the proof that the spins are the same. 
Let p{s)®1rxr (resp. p'(s) €5 l^'xi?') the orthogonal component of p{s) (resp. p'{s)) along 
ax(E)lRxR (resp. (Ta; "Xi 1 i?/ x i?/ ) • Using equation (|72|) . equality H7U|I with v = x can be written 
as : 

Tr {Pn, {t)p{.s)) R = Tr (Pjy; {t)p'{s)) R' . (74) 
Since we have assumed by contradiction R> R', from (|74ll we have for every t: 

Tr (Pjvi {t)p{s)) < Tr (Pjv; {t)p'{s)) . (75) 

Now we will derive a contradiction by evaluating the previous inequality at t = and i = tt 
and using (|73|l . In fact we have: 

Tr(F^,(0)p(s)) <rr(Pjv;(0)p'(s)), 

thus 

rr(Pjv,(7r)p(s)) = -TriP^Mpi^)) > -Tr (Pjvj (O)p'(s)) = Tr (Pjvj (7r)p'(s)) . 
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The previous inequality contradicts equation (|75|l . Thus we conclude that h = I'l, which 
implies that equation (|^ holds. 

Since the two equivalent models S and S' have the same spin, the positive constant a 
in (|69|l is equal to one. Therefore, for every pair d,k G {l,...,ri}, Jdk and J^^ only differ 
possibly by the a sign factor. Using the same argument as in the main Theorem of [2| we 
can in fact conclude that there are only two possible case: The case where Jdk = Jdk ^'^'^ 
every pair d,k and the case where Jdk — ~Jdk Gvery pair rf, fc. If we are in the first 
case, then from the observability (which follows from controllability) of the model, we must 
have po = Poi ^^^^ equation H59|) holds. This would be the case a of part 4. of Condition 
(*). On the other hand, if J^^ = —Jkd for every pair 1 < k < d < we may conclude 
using Theorem n In fact we consider the homomorphism (j) given by the Cartan involution 
associated to the odd-even decomposition as in the previous part of the proof. Conditions 
(PJ hold, thus, since the models are equivalent and observable, we get that: 

thus equation holds. This concludes the proof of the Theorem. 

6 Proofs of Theorems |4HS1 

In the proofs of Theorems ^ and |S1 we shall use the following two types of elementary k xk 
matrices: 

Cfe :=dm<?(-l,l,-l,...,(-l)'=), Tfc==adzag(l,l,l,...,l). (76) 

The matrix Ck is diagonal with alternating elements while Tk is antidiagonal with all ones 
on the secondary diagonal and zeros everywhere else. Obvious properties of these matrices 
are the following 

Cl = Tl = lky.k. Tk=T^. (77) 

We are interested in the action of these matrices by similarity transformation on diagonal 
and tridiagonal k x k matrices. In particular, let us denote by D a generic, real, diagonal, 
k X k matrix and by a generic, real, kxk, tridiagonal matrix, which is also symmetric and 
it has zero diagonal. If M° denotes the antitransposed of M, namely the matrix obtained 
by reflecting about the secondary diagonal, we can easily verify the following properties. 

1. 

CkDCk = -D, CkFCk = -F, (78) 

2. 

TkDTk = D\ TkFTk = F\ (79) 
Now we are ready to prove Theorems 0] and 

Proof of Theorem |4| The matrices icz and tax have (for every value of the spin) the 
following structure 

(81) 
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where F and D have the structure above specified with fc := and P is a A: x fc real matrix 
of all zeros except in the (fc, l)-th position. Now use 

which is orthogonal and therefore unitary. 

We calculate, using the first ones of H78|l and (|79|l 

ta, := Uia.U* = I (^^ J^^^ . (83) 
Moreover, using the second ones of (|78|l and (|79|l . we have 

la, -.^ Uia^U* ^ i (^^^f^^ (84) 

It is easily seen that iaz and iax are symplectic, by observing that CkPTk is a real symmetric 
matrix (only the (A;,A:)-th element is different from zero). Therefore sp{^) contains a 
subalgebra conjugate to the one spanned by ia^ and iaz and therefore iay and the Theorem 
is proved. □ 

We now proceed to the proof of Theorem 

Proof of Theorem |5l In this case we set k := ^^^^ ■ The matrix iaz has the form 

-L» 

iaz:^i\ | , (85) 
D" 

with D of dimension k x k. Moreover iax has the form 

F V 

iax ■■=i\ I , (86) 

w F" ^ 

where F is as above and v (w) is a vector of dimension k with only the last (the first) 
component different from zero, and the components different from zero are equal for v and 
w. We use the unitary matrix 

f/:= ( 1 ^0 , (87) 

which is easily seen to be unitary by H76|) (|77|l . We calculate. 

/ \{TkD-Tk - CuDCu) ^{{-ifTkD-Cu - CuDTk) 
iaz UiazU* = i\ 

\^{TkDCk^{-l)''CkD^Tk) ^{CkD-Ck-TkDTk) 

Using the first ones of H78() and (|79|l . we find that the diagonal blocks are zero. Moreover, 
the remaining elements of the matrix are real so that iaz is real. Analogously, we calculate 

iax ■■= UiaxU* = (89) 
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* * ^{TkFTk + CkF'^Ck) 

where we have denoted by a star * the components that can be obtained from the requirement 
that the matrix is skew Hermitian. Now the (1, 1) and (3, 3) blocks are zero from the second 
ones of properties (|75)l and l(7!?|) . while the (2,3) block is zero because of the structure of the 
vectors v and w. All the other blocks are purely real matrices so that idx is also in so{n) 
and this completes the proof. □ 

We now give the proof of the negative result in Theorem El 

Proof of Theorem |S| Assume that there exists a matrix X e SU{Nj) such that 

(90) 
(91) 
(92) 

with R^, Ry and R^ in so{Nj). Then we can use the AI Cartan decomposition of SU{Nj) 
[5] to write X as 

X ^ K1AK2, (93) 
with Ki and K2 in SO{Nj) and A diagonal i.e. 

A:=diag{e"l'\...,e ^). (94) 

Therefore we can write 

KiAK2t<J,,y^,K^AKl = Rx^y,,, (95) 
or, defining Rx.y,z '■— R^^y^zKi which is also real skew-symmetric, we can write 



XlGxX* 


'■= Rx, 


XiayX* 




Xia^X* 


■■^Rz, 



K2i<Jx,v,zKl = ARx^y^.A. (96) 



The real matrices Rx,y,z must satisfy the same basic commutation relations (|35f) of 



iay and iaz and have the same eigenvalues of iax, i<Jy and iuz, namely for a (half integer) 
spin j, ±ji, ±(j + l)i,...,±ii. We now study the structure of Rx,y,z in H95|l and get a 
contradiction with these facts. 

First notice that, since A is diagonal as in (|94|l . the action of A on the right hand side 
of H96|) namely, R — > ARA changes the (real) clement Vjk of R into rjke~'^^'^^~'^''\ Since 
the entries on the left hand side of l|96|) are cither all purely imaginary or all purely real, if 

— (t>k is not a multiple of then we must have rjk = 0. Consider the indices 1, Nj and 
let O be the set of indices k such that 0i — 0^ is and odd multiple of ^ and £ the set of 
indices k such that — (j)k is an even multiple of ^ , and Af the set of indices k such that 

— 0fc is not an integer multiple of ^ 

From H9t)|) it follows that since ia-y is real the terms rjk of Ry where j and k belong to 
different sets must be zero because in that case e^^'^^~'^''^ in H94|l has a nonzero imaginary 
part in this case. Therefore only the elements rjk where j and k belongs both to O, or both 
to £, or both to Af, are possibly different from zero. Therefore after possibly a reordering 
of rows and columns, which corresponds to a similarity transformation by a permutation 
matrix, Ry must be of block diagonal form and without loss of generality and for simplicity 
we shall assume only two blocks (rather than three). Therefore we write 



22 
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where Yn has dimensions Uo x Hq with Uo the cardinahty of O and Y22 has dimension 
(n — rio) X (n — rio). Both Yn and I22 are skew-symmetric matrices. An analogous argument 
shows that, after possibly the same re-ordering of column and row indices, Rz can be written 
as 

where Z12 is a general matrix of dimensions rio x (n — no). 

Now consider the possible values for Uo- Uo odd is to be excluded because this would 
cause detVii = in H97|l and this contradicts the fact that Ry has no zero eigenvalues. 
Moreover no ^ (n — Uo) (i.e. rio ^ ^) would cause R^ to have determinant equal to zero. 
This can be easily verified by calculating 

det{Rl) = {detRzf = det ^^^^^^^ -ztzi2 ) ' ^^^^ 

since, in this case, at least one of the matrices on the diagonal blocks does not have full rank. 
These considerations already exclude the cases where is an odd number as for spins ^, 

|, I etc. and we can assume Ry and Rz of the form l|?l7jl and with rio = To obtain 
a contradiction in this case too we first notice that, since Yn and have even dimension 

and are skew-symmetric, we can apply a similarity transformation T -.^ [ ^ ) j with 



Ti and T2 orthogonal so that TRyT^ is block diagonal 



TRyT^ ^{Di,D2,...,Dn^), (100) 

2 

where the 2x2 block Dk has the form 

^-=(-1 o)' 

where each Ik corresponds to a couple of complex conjugate eigenvalues of Ry so that Ik — 2 
with p odd corresponds to the pair ±|i. Moreover we choose T so that the first ^ blocks 
are ordered according to the increasing value of Ik and the same thing for the last blocks. 
We shall therefore assume this structure of Ry in the remainder of the proof. We notice 
also that the transformation TRzT^ does not change the structure of Rz as Z12 in 1)98(1 
was chosen to be a general ^ x -j- real matrix. Express Z12 in terms of 2 x 2 blocks L/fc, 
f,k — 1, fc = + 1, which is possible since ^ is an even number. Now, we 
impose the fact that Ry and Rz have to satisfy the same commutation relations as i<Ty and 
iaz- In particular, we must have 

[[Ry,Rz],Ry]^ Rz- (102) 
This equation gives for the Lfk block the following one 

2DfLfkDk - LfkDl - DjLfk = Lfk. (103) 



If we write the generic L fk as 



and recall the structure oi Dj and Dk, 
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we obtain the following equations for 02 and 03 (and analogous equations for ai and 04) 



2lklfa3 = {l-l}^ll)a2, (106) 

2lklfa2 = il-l}~ll)a3. (107) 

Combining these, we obtain 

4lll}as^{l~ll~l}fa,, (108) 

which shows, taking the square root of both sides, that the only possibilities to have 03 and 
therefore 02 different from zero are the cases + Zfc = ±1. In these cases we can easily see 
that 

as = -02. (109) 

Similarly, one finds that we have 04 and oi in (|104|1 different from zero if and only if 
If + Ik = il and, in these cases we have 

ai — 04. (110) 

In conclusion, all the blocks Lfk are zero except the ones corresponding to indices / and k 
with neighboring values of // and Ik which have the structure 

^/^-(-; i)- (111) 

Therefore Rz has the form in (|98|l where the /— th block row of Z12 has at most two blocks 
different from zero and with the structure in Hlll() . We denote these blocks by Pf and 
Sf, where P (S) stands for 'predecessor' ('successor') and correspond to the index k such 
that Ik — If — 1 and Ik ~ If + 1, respectively. Now, we argue that a matrix with this 
structure must necessarily have all the (purely imaginary) eigenvalues with multiplicity at 
least two and this gives the desired contradiction because Rz should have the same spectrum 
of i<Jx,y,z which consists of all simple eigenvalues. In order to see this fact, re-consider the 
block structure of Ry in (|100|l . If the blocks corresponding to eigenvalues ±ii and ±|z 
belong to the same half, then the corresponding matrix Rz will have a two dimensional 
block row (or column) equal to zero, and therefore will be an eigenvalue with multiplicity 
at least 2. Therefore we can assume that these two blocks belong to two different halves 
and by the ordering we have imposed they must be the first ones of each half. Assume 
that the block corresponding to ±ij is in the first half. If this is not the case, consider 
the transposed of Rz and repeat the arguments that follow. It is possible to choose a block 
diagonal similarity transformation 

C/:=d^a5(Gl,G2,...,G«^,Fl,F2,...,F«^), (112) 

4 4 

with all the G/'s and F/'s 2 x 2 orthogonal matrices so that URzU'^ has the same structure 
as before but all the matrices Pj and Sj are scalar matrices. We construct the matrix U 
proceeding by block rows. The first block row contains only Si as ^ has no predecessors. 
All the zero blocks remain zero and Si is transformed as 

GiSiF^. (113) 

We choose Fi = l2x2 and Gi, which has the general form 

V —sin[u) cos[o) J 
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so that sin{9)x + cos(6)y = if 5*1 := I ^ ^ | . This will give a scalar matrix. At the 

\-y X J 

generic /— th block row, we have, at the most two nonzero blocks Pfk and Sfb where we 
now use an extra index k and b to indicate the block column to which they belong. They 
transform as 

Pfk ^ G/F/fcFj, (115) 

and 

^fb ^ GfSfbF^ , (116) 

respectively, while all the other blocks remain zero. If Fk has not been chosen before we 
set Fk = l2x2- In any case, we choose G/ as before to make GfPfkF^ a scalar matrix. 
We then choose Fi, to make GfSfbF^ a scalar matrix. Gf and Ff, had not been chosen at 
previous steps. This is obvious for Gf and follows by an induction argument for F^ since 
all the F matrices chosen before the /— th step correspond to predecessors and successors 
with (column) indices strictly less then h (recall that in the two halves of the matrix Ry 
the blocks are arranged in increasing order of (absolute value of) eigenvalue). In conclusion, 
modulo the similarity transformation defined by U in 11121) , we can assume that has the 
form 

Rz=K®h^2. (117) 

where K \s a. skew-symmetric x ^ matrix. By known results on the eigenvalues of the 
Kronecker products of two matrices, it follows that the eigenvalues of Rz are the same as 
those of K each with multiplicity at least two. This gives the desired contradiction and 
concludes the proof of the Theorem. □ 



7 Conclusions 

This paper has presented a collection of mathematical results concerning the input-output 
equivalence of quantum systems. Models that are equivalent cannot be distinguished by an 
external observer and therefore the determination of parameters in a quantum Hamiltonian 
can only be obtained up to equivalent models. Motivated by recent results on the isospectral- 
ity of quantum Hamiltonians |14| in molecular magnets, we have completely characterized 
the classes of spin networks which are equivalent. In several cases, the characterization of 
equivalent models can be obtained through a Lie algebra homomorphism which is suggested 
by a Cartan structure of the underlying dynamics. 

We believe many of the results and the concepts presented in this paper for quantum 
systems could be generalized to classes of systems relevant in other applications with both 
dynamics and output linear in the state. This will be the subject of further research. 
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Appendix A: Additional results and proofs 

The proof of the following Lemma can be obtained with a formal modification of the proof 
of Lemma 4.4 in |2| and it is therefore omitted. 

Lemma 7.1 Let S and S' be two equivalent models models. If W and W' are two given 
Hermitian matrices such that 

Tr{Wp{t)) = Tr{W'p'{t)), (118) 
for every pair of corresponding trajectories p{t) and p'{t), then it also holds 

Tr([W,lMt))=Tr{[W',l[^)^]p\t)), V fc G {1, . . . , n}, V« £ {a:, 2}, (119) 
up to a permutation of the indices ^. 

''This permutation is the same and fixed for all the results where it is mentioned 
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Proof of Lemma 15.11 

We first state a Lemma whose proof can be obtained from the proof of Lemma 5.1 in |2] 
and then proceed to the proof of Lemma lOI 

Lemma 7.2 Assume that for all t > 0, all the possible trajectories p{t) of S and corre- 
sponding p'{t) of S', for fixed values 1 < ki, . . . ,kr < n, and fixed Vj G {x, y, z}, we have: 



Tr(/fci„i,...,fc,^,p(i)) =Tr (^l'k,vi,...,k,.v^P' (t)) ■ 

Then: 

1. equation (|120|l holds for any possible choice of the values of £ {x,y,z}; 
2. 



(120) 



il'j , [iK , A' 

dvj' I- kvT. 



; lkiV\ ,...,kj.V 
' ^ki vi ,. . . ,krVj~ 



p{t) 
p'{t) 



(121) 



for all the indices 1 < d ^ k < n and every {wj ^ v^} € {x, y, z}. 



We now proceed to the proof of Lemma 15.11 First notice that from Lemma 17.21 it is 
enough to prove Ht)5|) and (|66|l for a particular choice of {vj} and v. Moreover, we have, for 
d> k, 

[^hz, A]] - -Jkd^hz,dx- (122) 

1. By applying Lemma [^| (equation (|121|l ) to H64|l and using H122|) we get: 

PTr{[-JJ,^iIJ,^^^^,h,v,^...MrV^,p{t)) = (iTr(^ -J-kd^hzM^h^v^^.^k^v, p'{t)) ■ (123) 
We may assume, without loss of generality, that k — kj and Vj = x. In this case we have: 

^Jkd \_-^kz,dxT ^kivi,...,krv,^^ = Jkd'HkiVi ,. . . ,kjy ,...,krV,. ,dx- 

Combining the previous equality with (|123|l . equation H65|) follows easily. 

2. Using the same procedure, we obtain again equation 11231) . but now both indices k and 
d are in {fci, . . . , k^}. Assume, for example that ki — k and k2 = d, and take Wfcj = = x. 

Now we have: 

[Ikiz,k2X7 Ikix,k2X,...,krV,.] — -^kiy ,k2X^ ,k3V3 ...,k,.Vr ^ (124) 

where, with this notation, we mean that in the position we have the matrix ct^. Thus, 
combining equations (|123ll and 112411 . we get: 



PJkik2Tr {lkiyMx^Mv3...,krV^P{t)) — /^'^/fcifca^'^ (^hiyMx^. 

Using the same procedure, we conclude: 

fiJkik2Tr {lkiy,k2y^Mv3...,k^v^p{t)) = /3'4ife2^^ {^lyM 



k3V3...,krV., 



.P'{t) 



y^ ,k3V3...,krVr 



.Jit) 



and 



(125) 
(126) 
(127) 



PJktk2Tr {lkiy,k2Z^.k3V3....k^V^Pit)) = /J'Jfcifc^Tr [lk^y^k2Z^ Mv3-.k^VrP' 

Adding together equations p25|l . (|126|l . and (|127|) and using ((S^ . we get: 

(i{lk2{lk2 + l))Jk,k2Tr[h,y,k3V3-,KvAt)) = /?'(?;, (4, + l))4,fc,Tr [h,y,k3V3...,k^vJ{t)) , 

as desired. 
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Proof of Lemma 15.31 

Proof. 

We recall the formulas 172|) and (|73|l to be proved, i.e.: 

(e'"-* ® lflxfl)/i.(e-'"^* ® Iflxfl) := PjVi W li?.xii, (128) 

where the matrix P(-) is periodic with period 27r, and 

Fw,(^) = -Pjvi(0) = a.. (129) 

The proof can be done directly by computing the matrix above. This is simplified by the 
fact that the matrix is always a diagonal matrix. We will give an outline of the argument 
when li is half integer spin. The idea is to use the representations for the matrices az and ax 
given by equations H8U|I and H81|) . The case of integer spin can be derived similarly starting 
with the representations given by equations H85|l and (|86|l . 
Using equations and (jH^I, we obtain 



e 



The properties of the matrices D, P and F are described in Section |B1 Moreover D — 
diag{j,j — 1, ^) for a half integer spin j. By using these properties, it follows that all the 
time depending terms in equation H13U|) are of the form e**. Thus matrix (|130|l is periodic of 
period 27r. The fact that the dependence is of the type e**, in turn, implies that equations 
ijnSj) and (inni hold. 

□ 
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